Molecular modeling and simulation approaches to characterize potential molecular targets for burdock inulin to instigate protection against autoimmune diseases

In the current study, we utilized molecular modeling and simulation approaches to define putative potential molecular targets for Burdock Inulin, including inflammatory proteins such as iNOS, COX-2, TNF-alpha, IL-6, and IL-1β. Molecular docking results revealed potential interactions and good binding affinity for these targets; however, IL-1β, COX-2, and iNOS were identified as the best targets for Inulin. Molecular simulation-based stability assessment demonstrated that inulin could primarily target iNOS and may also supplementarily target COX-2 and IL-1β during DSS-induced colitis to reduce the role of these inflammatory mechanisms. Furthermore, residual flexibility, hydrogen bonding, and structural packing were reported with uniform trajectories, showing no significant perturbation throughout the simulation. The protein motions within the simulation trajectories were clustered using principal component analysis (PCA). The IL-1β–Inulin complex, approximately 70% of the total motion was attributed to the first three eigenvectors, while the remaining motion was contributed by the remaining eigenvectors. In contrast, for the COX2–Inulin complex, 75% of the total motion was attributed to the eigenvectors. Furthermore, in the iNOS–Inulin complex, the first three eigenvectors contributed to 60% of the total motion. Furthermore, the iNOS–Inulin complex contributed 60% to the total motion through the first three eigenvectors. To explore thermodynamically favorable changes upon mutation, motion mode analysis was carried out. The Free Energy Landscape (FEL) results demonstrated that the IL-1β–Inulin achieved a single conformation with the lowest energy, while COX2–Inulin and iNOS–Inulin exhibited two lowest-energy conformations each. IL-1β–Inulin and COX2–Inulin displayed total binding free energies of − 27.76 kcal/mol and − 37.78 kcal/mol, respectively, while iNOS–Inulin demonstrated the best binding free energy results at − 45.89 kcal/mol. This indicates a stronger pharmacological potential of iNOS than the other two complexes. Thus, further experiments are needed to use inulin to target iNOS and reduce DSS-induced colitis and other autoimmune diseases.


Preparation of ligands and receptors
Initially, we downloaded the crystal structures of iNOS (PDB ID; 3HR4), COX-2 (PDB ID; 5IKQ), TNF-alpha (PDB ID; 2AZ5), IL6 (PDB ID; 7NXZ), and IL-1β (PDB ID; 4G6J) from the RCS PDB database.Subsequently, these structures were processed for further analysis by removing any extra attached water molecules and were saved in .pdbqtformat 26 .The ligand structures were obtained from the PubChem database, and their energy was minimized.Subsequently, the resulting structure was subjected to an energy minimization algorithm using the Swiss-PdbViewer v.4.10 software program 27 .To ensure the protein was suitable for docking with AutoDock Vina, we first obtained its three-dimensional structure and meticulously refined it.This involved addressing any missing atoms or loops within the protein.We then optimized the protein structure by eliminating water molecules and introducing essential hydrogen atoms.The refined protein structure was saved in the PDBQT file format, which www.nature.com/scientificreports/AutoDock Vina specifically requires for protein input.Similarly, we acquired the three-dimensional structure of the ligand and followed a comparable procedure for its preparation.The ligand underwent optimization by adding hydrogen atoms and eliminating unfavorable conformations.Additionally, partial charges were assigned to the ligand atoms, commonly employing force field-based methods.Finally, the prepared ligand structure was saved as a PDBQT file format, ready to be docked with the protein using AutoDock Vina.

Molecular docking
For protein-ligand docking, various computational strategies are routinely employed, especially for discovering small drug-like molecules.In this study, we utilized AutoDock Vina to examine the binding interactions between Inulin and iNOS, COX-2, TNF-alpha, IL6, and IL-1β 28 .AutoDock Vina employs an exhaustiveness parameter, which determines the thoroughness of the search algorithm during the docking process.The amino acid residues that constitute the active site of the protein were firstly delineated, which a grid box of size x = 15.1651137855Å, y = 19.301165182Å, and z = 13.9173764782Å, with a center dimension of x = 18.4802485932, y = -18.9999491012,and z = -50.2219055were set to define the active site.Conversely, a blind docking protocol was employed to study the interactions, the grid box with size x = 62.2440970136Å, y = 77.8606218435Å, and z = 62.064913585Å, with a center dimension of x = 31.0115192852,y = 71.8214808942,and z = 80.3147800551 was found to cover every residue of the protein.The docking procedure was then executed using an exhaustiveness of 8.This allows for a more exhaustive exploration of the binding landscape and increases the chances of identifying optimal binding modes.In our approach, we set the exhaustiveness parameter to 16 to enhance binding accuracy.For each case, 10 different poses were initially generated, and based on the best docking score, the best conformation was selected to examine the protein-ligand binding mechanism using PyMol [29][30][31] .

Molecular dynamics simulation
To understand the structural inhibition effect of Inulin on iNOS, COX-2, and IL-1β cytokines, the docked complexes were inspected by a molecular dynamics' simulation procedure using Amber 20 simulation software.We initially employed the "FF14SB" force field to define the protein 32 , and then added a cubic box of the TIP3P water model (with a buffer distance of 10 Å) to solvate the protein.Both negative and positive counter ions (Na + and Cl − ) were further added to balance the system 33 .To ensure the system is fit for optimal performance, a total of 1000 energy minimization steps were performed including the 500 steepest descent and 500 conjugate gradient steps.The system was steadily heated from 0 to 325 (K) 34 and the hydrogen bonds are subsequently controlled via the SHAKE method.During the MD simulations procedure, the Langevin dynamics approach was used i.e., with pressure of 1 atm and temperature of 310 K 35 .At a constant pressure and temperature, each system was first equilibrated for 10 ns (ns) time.Following that, the MD simulation was conducted for 150 ns (ns) total time at a constant temperature of 325 K 36 .Finally, data visualization and analysis were performed with the help of Origin Software.We further thoroughly examined the conformational alterations that eventuated during simulations.To analyze these alterations, we utilized CPPTRAJ module within AMBER20 to investigate several factors namely root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), Radius of gyration (RoG) and Principal component analysis (PCA).The Radius of gyration (RoG) is a valuable measure for analyzing protein size and compaction 37 .

Binding free energy calculation
To evaluate the binding strength in molecular simulation, we employed the MM/PBSA.pyscript to calculate the binding free energy 38,39 .This method, known for its broad applicability, offered several advantages over other approaches, including computational efficiency and reduced time requirements.For each complex, the van der Waals (vdW), electrostatic, solvent-accessible surface area (SASA), and as well as the generalized Born (G.B.) components were determined on the basis of the entire simulation trajectory.This approach combines molecular mechanics simulations, which describe the interactions between atoms, with implicit solvent models, which describe the interactions between the protein and solvent [40][41][42][43][44] .Hence, we also applied this approach here to accurately compute the total binding free energy of the protein-ligand complexes.Mathematically the binding free energy can be estimated as: Different contributing components of total binding energy were calculated by the following equation: It has a wide range of applications i.e., used to estimate the binding energy for proteins in different studies including SARS-CoV-2 and neurological disorders 43,[45][46][47][48][49] .

Principal component analysis (PCA)
Each clustering of the trajectory was created on the basis of calculation and diagonalization of the covariant atomic fluctuation using the principal component analysis (PCA) 50,51 .In PCA analysis both diagonal eigenvectors and analogous eigenvalues values were extracted.The principal components (PCs) that provide the information on atom motions and displacement are also known as eigenvectors and has already been discussed previously 52,53 .The following equation was used for the cluster trajectory analysis. (1) X denotes the response of two principal components, kB is Boltzmann constant, while the P(X) represents the dispersion of the framework's likelihood.The Free energy landscape (FEL) was built from the two PCs.

Results and discussion
In the current study, molecular modeling and all-atoms simulation approaches have been integrated to computationally validate the potential target for Inulin, which is an active ingredient of Burdock.Reduced expression of IL-6, IL-1β, TNF-alpha, iNOS and COX-2 in the experimental analysis of burdock extract revealed these receptors as the potential targets in immunocompromised diseases and particularly DSS-induced colitis [54][55][56][57] .We employed structural and simulation methods to investigate the binding of Inulin with these targets and further validated by using binding free energy approach.The docking results demonstrated that the iNOS-Inulin exhibit the best docking energy of − 7.54 kcal/mol followed by IL-1β with the docking score of − 5.08 kcal/mol.Moreover, the docking results for COX-2 and TNF-alpha complexes revealed to be − 4.80 kcal/mol, respectively.For IL-6, the docking energy was calculated to be − 4.00 kcal/mol, which is the lowest among all.The biomolecular interaction paradigm of each complex was explored and the interaction features were comprehensively explored.For instance, the iNOS-Inulin complex reported various group of interactions including van der Waals, hydrogen bonds, hydrophobic and charge interactions.The iNOS-Inulin complex stands the best with 8 hydrogen bonds established with Trp194, Arg199, Ile201, Gly202, Ser242, Phe369, Asn370, Gly371 and Trp372, while the salt-bridges, pie-anion, alkyl, pie-alkyl and attractive charges contacts were also frequently observed.The 3D binding mode of inulin the binding cavity of iNOS is given in Fig. 1a, while the 2D detailed interaction pattern of inulin with iNOS is shown in Fig. 1b.The second-best docking complex was IL-1B-Inulin, which established two hydrogen bonds and multiple hydrophobic interactions.Among the hydrogen bonds Glu25 and Lys74 residues while the Lys74, Lys77 and Pro131 residues were involved in hydrophobic interactions.The residues Lys74 and Lys77 were particularly involved in salt bridges while Pro131 connected with inulin through an alkyl interaction.For instance, these residues are also reported to act as inhibitor hotspots for other drug targets 58 .The 3D binding mode of inulin the binding cavity of IL-1B is given in Fig. 2a, while the 2D detailed interaction pattern of inulin with IL-1B is shown in Fig. 2b.
Moreover, the COX2 and TNF-alpha demonstrated a similar docking score of − 4.80 kcal/mol.The interaction paradigm of COX2-Inulin revealed a single hydrogen bond with Asn375, while the pie-anion interaction was established by Phr209 and alkyl interaction was reported with Val228 and Val349.Although the complex demonstrated relatively good docking score but the interactions of the ligand with the active site residues are relatively lower in numbers.Interestingly our results strongly corroborate with previous literatures where blocking of these residues in the COX2 and TNF receptor led to the inhibition of the pathological role of these proteins [59][60][61] .The 3D binding mode of inulin the binding cavity of COX2 is given in Fig. 3a, while the 2D detailed interaction pattern of inulin with COX2 is shown in Fig. 3b.
On the other hand, the TNF-alpha bound with inulin reported two stronger hydrogen bonds with Tyr119 and Gln149, while three salt bridges by Tyr59 were observed.Although the number of bonds is more in number than the previous compounds but still the binding energy align with each other.The 3D binding mode of inulin the binding cavity of TNF-alpha is given in Fig. 4a, while the 2D detailed interaction pattern of inulin with TNF-alpha is shown in Fig. 4b.
To see how the binding of inulin occur with IL-6 the docking was also performed against the active site residues of IL-6.The docking results revealed that Leu1 and Lys185 established robust hydrogen bonds with inulin, while the salt bridges were reported with Glu97 and Lys182 residues.Lys182 was also observed to have a hydrogen bond with the backbone of inulin, while the tails are attached by Leu1 and Lys185.The backbone of the inulin is strongly attracted inside the cavity by the salt bridges.In many other diseases such as COVID-19 IL-6 serve as a target to rescue the host immune response where direct inhibition of these demonstrated the reduced role of IL-6 in viral pathogenesis 62 .Since, our results also block the same residues thus may induce more potent pharmacological properties.The 3D binding mode of inulin the binding cavity of IL-6 is given in Fig. 5a while the 2D detailed interaction pattern of inulin with IL-6 is shown in Fig. 5b.

Root mean square deviation
Assessing the dynamic stability of a protein bound to a ligand is essential to demonstrate the pharmacological activity of that particular compound.The compound IL-1β-Inulin stabilized at 1.3 Å.In the case of compound   IL-1β-Inulin the RMSD initially remained stable with no major deviation.No significant deviation was observed during the simulation except a minor increment in the RMSD graph at 0-10 and 85-90 ns.Then the RMSD gradually decreased and increased following this pattern until the end of 150 ns simulation.However, the complex demonstrated stable dynamics with an average RMSD of 1.30 Å, because when RMSD ≤ 2 Å can prove the reliability of the docking method.For the COX2-Inulin, the RMSD increased gradually since the beginning.The RMSD had large fluctuations and continues to increase with abrupt increase and decrease between 30 and 60 ns.The COX2-Inulin reported significant structural perturbation after 30 ns, the maximum RMSD value is 4.8 Å.The complex demonstrated a uniform RMSD pattern between 60 and 100 ns with the gradual increasing trend but then significant perturbations were seen in the complex between 100 and 105 ns, the maximum RMSD value is 4.2 Å.Then maintain a steady trend.Because the RMSD fluctuates greatly, the difference between the maximum value and the minimum value is close to 4.0 Å, and the complex demonstrated unstable dynamics with an average RMSD of 3.5 Å.When RMSD ≥ 2 Å can't prove the reliability of the docking method.In general, the compound IL-1β-Inulin is relatively stable, but COX2-Inulin is unstable.On the other hand, the iNOS-Inulin complex that reported the best docking score among the all with highest number of hydrogen bonds and other interactions also demonstrated stabilized dynamic behavior throughout the simulation.The complex stabilized at 3.0 Å and continues to remain uniform with minor deviation at 75 ns.The complex overall attained a very stable dynamic with an average RMSD of 2.8 Å until the end of simulation.The RMSD for the iNOS, COX-2, and IL-1β Apo forms is shown in Fig. 6.Dynamic stability of the ligand-bound complexes always yields better pharmacological output as witnessed by many studies, and hence our results strongly align with the previous results were stable complexes of the aforementioned results reported stronger inhibition than the less stable 45,46,58,62 .Overall, this shows that inulin could primarily target iNOS and may supplementarily target the COX-2 and IL-1β during  www.nature.com/scientificreports/ the DSS-induced colitis to reduce the role of these inflammatory machines.The RMSD results for the top three complexes are given in Fig. 6a-c.

Residual flexibility and hydrogen bond analysis
Root Mean Square Fluctuation (RMSF) analysis was performed to investigate the dynamic behavior of the protein structure throughout the simulation trajectory.The RMSF values were calculated for each residue in the protein.Figure 7a shows the RMSF profile of the protein over specific residues.The x-axis represents the residue number, while the y-axis represents the RMSF values in angstroms.The RMSF values were calculated by measuring the deviation of each residue's position from its average position during the simulation.The RMSF analysis of IL-1β-Inulin revealed distinct regions of flexibility within the protein structure.In the case of iNOS-Inulin complex, residues 80-100 demonstrated higher fluctuations while the rest of the regions reported minimal fluctuations.
Hydrogen bonding calculation is one of the key assessments that help in determining the pharmacological potential of a drug/inhibitor.It is an essential approach to reveal the potency and binding strength of the interacting molecules.This approach has been widely applied to understand the pharmacological mechanism of a particular drug, and the interaction mechanism of two or more proteins to reveal the mechanism of a disease or bio-catalytic process [63][64][65][66][67] .Then we performed hydrogen bond analysis to investigate the dynamic behavior and intermolecular interactions within the protein structure throughout the simulation trajectory.Hydrogen bonds were calculated based on distance and angle criteria using a cutoff of 3.5 Å for the donor-acceptor distance and 30° for the donor-hydrogen-acceptor angle.The analysis revealed that all the complexes (IL-1β-Inulin, COX2-Inulin, iNOS-Inulin) maintained several hydrogen bonds within the protein structure in association with the ligand Fig. 7b.Each line represents a hydrogen bond, with the donor and acceptor residues labeled accordingly.Interestingly, in case of IL-1β-Inulin residues involved in hydrogen bonding interactions were found to be consistent with the regions of high flexibility identified through the RMSF analysis.Residues 20-25, 42-53, 62-70, and 100-130, which exhibited higher RMSF values, also showed a higher propensity for forming hydrogen bonds with neighboring residues.

Radius of gyration
Protein packing determination or protein size estimation demonstrates information regarding the essential events that occurred during the molecular simulation.We estimated the structural compactness of each complex during simulation by calculating the radius of gyration (Rg) as a function of time.At the beginning, the Rg of the compound IL-1β-Inulin remained stable, the Rg continued to increase with abrupt increase and decrease between 20 and 30 ns, and there was a maximum value between 25 and 30 ns, with a maximum value of 15.3 Å.The Rg remained stable with no major deviation between 30 and 65 ns.Starting from 65 ns, the Rg pattern also reported a gradual decrease in the Rg trajectory.The Rg started from 15.00 Å, and reached to a maximum of 14.80 Å at 150 ns.An average Rg was calculated to 15.00 Å, thus convey a stable binding of the ligand with no significantly increase or decrease in the protein size during simulation.And the Rg image of the COX2-Inulin is not very stable, showing significant fluctuations overall.The Rg for the COX2-Inulin complex started to decrease initially and reached 19.20 Å at 5 ns and then equilibrated.A straight uniform Rg was seen until 35 ns and then an abrupt increase was observed in the Rg pattern, and reached to a maximum of 19.8 Å. Subsequent www.nature.com/scientificreports/ Rg also fluctuated at 60 ns, 100 ns, and 125 ns.This shows significant structural perturbation occurred during the simulation and thus causes binding and unbinding of the ligand in the pocket, and an average Rg was calculated to 19.30 Å.Moreover, the Rg for iNOS-Inulin started to increase initially and reached to 18.8 Å, however they started to decrease gradually and attained the tighter structural packing at 30 ns.The Rg continues to follow the same pattern until 130 ns but a small perturbation was observed for a shorter period and then again decreased back.An average Rg for the iNOS-Inulin complex was calculated to be 18.45 Å.The Rg results for the top complexes are given in Fig. 8a-c.

Clustering of protein's motion
The protein motions within the simulation trajectories were clustered utilizing principal component analysis (PCA).Each system structural modifications brought about by mutation were plotted using principal component analysis (PCA).Utilizing a trajectory of 0-150 ns MD simulation, PCA was performed to learn about the conformational states of the IL-1β-Inulin, COX2-Inulin and iNOS-Inulin complexes.The eigenvalues of the covariance matrix served as the foundation for conducting principal component analysis (PCA), enabling the assessment of the overall motion of Cα-atoms in the hits and complexes.The motions observed in two dimensions were represented by PC1 and PC2. Figure 9a-c illustrates the distributed principal components for each complex.The dots represent the various frames, which progress from orange to indigo.The IL-1β-Inulin approximately 70.0% of the total motion was attributed to the first three eigenvectors, while the remaining motion was contributed by the remaining eigenvectors.In contrast, for the COX2-Inulin eigenvectors contributed 75.0% to the total motion.Furthermore, the iNOS-Inulin contributed 60.0% of the total motion by the first three eigenvectors.
The conformational transition between them can be clearly distinguished by the orange and indigo colors.As a result, the overall internal motion is reduced, leading to effective binding effects.
In order to explore thermodynamically favorable changes upon the mutation, motion mode analysis was carried out.PCA was applied using the covariance matrix of the 3D positional coordinates of different IL-1β-Inulin and COX2-Inulin complexes obtained from a 150 ns MD simulation.This allowed identification of the principal components (PCs) representing specific motion modes performed by these complex systems.The first PC captured the dominant motion observed in the analysis.Figure 9d-f demonstrate that the IL-1β-Inulin achieved a single conformation with the lowest energy, COX2-Inulin and iNOS-Inulin exhibited two lowest-energy conformations.

MM/GBSA analysis
Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) analysis was performed to assess the binding free energy and contribute to our understanding of the protein-ligand interactions.The MM/GBSA calculations were carried out for the protein-ligand complexes sampled during the molecular dynamic's simulations.Table 1 summarizes the calculated MM/GBSA energy components for each protein-ligand complex.The binding free energy (ΔG_bind) is represented as the sum of the contributions from different terms, including the molecular mechanics energy (ΔG_MM), the solvation energy (ΔG_GBSA), and the entropy contribution (ΔS).values suggest favorable binding, while positive values indicate unfavorable binding.In our study, IL-1β-Inulin and COX2-Inulin displayed ΔG bind values of − 27.76 kcal/mol and − 37.78 kcal/mol, while the iNOS-Inulin demonstrated the best binding free energy results, − 45.89 kcal/mol, in contrast to the first two complex, respectively, indicating strong binding affinities.Further, their individual energy contribution was calculated which shows the residues Trp194, Arg199, Ile201, Gly202, Ser242, Phe369, Asn370, Gly371 and Trp372.The total energy contribution of these residues showed correlation of − 011, − 0.016, − 0.17, − 0.047, − 0.85, − 0.072, − 0.48, − 0.44, and − 0.034 respectively.Overall, the MM/GBSA analysis provided valuable insights into the binding free energy and the key molecular interactions between the protein and ligand.The results provide significant insights into the inhibitory potential of these targets by BFO and are showing better results than others 7,58,68 .These findings enhance our understanding of the structure-activity relationships and guide future drug design efforts targeting the protein of interest.

Discussion
Pro-inflammatory cytokines exert a central regulatory role in the pathogenesis and advancement of DSS-induced colitis 69 .In our experimental study, we observed a significant increase in the mRNA expression levels of proinflammatory cytokines, including IL-1β, IL-6, and TNF-α, following DSS administration, indicative of colitis induction in mice.However, intriguingly, administration of inulin led to a remarkable reduction in the upregulation of these cytokines.This suggests that inulin exerts a suppressive effect on the inflammatory response associated with DSS-induced colitis.
Similarly, the expression of inducible nitric oxide synthase (iNOS) and cyclooxygenase-2 (Cox-2) markedly increased following administration of dextran sulfate sodium (DSS).However, varied concentrations of inulin exhibited a dose-dependent reduction in the expression of these proinflammatory cytokines.This suggests a  www.nature.com/scientificreports/significant mitigating effect of inulin on the inflammatory response induced by DSS treatment 25 .iNOS and COX-2 play crucial roles in inflammatory signaling pathways, triggering the synthesis of nitric oxide (NO) and prostaglandin E2 (PGE2) in the mucosa.This activation can contribute to an increase in inflammatory cytokines such as TNF-α, IL-1β, and IL-6, thereby intensifying inflammation and causing tissue damage 70,71 .
In a previous study, buddlejasaponin IV demonstrated the capacity to diminish the expression levels of proinflammatory proteins such as iNOS and COX-2, alongside the mRNA expression of cytokines TNF-α, IL-1β, and IL-6 in RAW264.7 cells.This effect was achieved through the inhibition of the NF-κB signaling pathway, elucidating its potential as a regulator of inflammation at the molecular level.In cellular processes, the transcription factor NF-κB is intricately regulated by its association with an inhibitory subunit, IκB, predominantly located in the cytoplasm in an inert state.Upon stimulation, IκB undergoes phosphorylation, triggering its degradation and facilitating the translocation of NF-κB into the nucleus.Once in the nucleus, NF-κB orchestrates the transcriptional activation of genes responsible for various cellular functions.Notably, the compound buddlejasaponin IV exhibited a concentration-dependent inhibition of NF-κB translocation, thereby modulating its activity in the nucleus 72 .In our investigation we examined the possibility that inulin inhibits NF-κb activity in a similar pattern.
Liu et al. conducted an investigation into the synergistic impact of inulin in combination with Lactobacillus rhamnosus using a murine model of DSS-induced colitis.Their study revealed the potent anti-inflammatory attributes of inulin.Elevated levels of proinflammatory cytokines, including IL-1β, IL-6, and TNF-α, were observed in the DSS-treated group, a response mitigated upon the administration of inulin combined with Lactobacillus rhamnosus.Moreover, this combination led to an enhancement in the abundance and diversity of intestinal microbiota.Notably, alterations in the microbiota composition may indirectly influence iNOS levels, indicating a potential mechanism underlying the observed anti-inflammatory effects.Inulin is fermented by gut bacteria in the colon, generating short-chain fatty acids (SCFAs), including metabolites.SCFAs, particularly, serve as energy sources for intestinal cells and regulate immune cell differentiation 73,74 .Additionally, in stressrecurrent inflammatory bowel disease (IBD), inulin was found to improve symptoms by modulating microbiota composition and reducing inflammation, which may involve the regulation of iNOS.Inulin decreases the serum inflammatory markers IL-6, ACALP and inflammatory cytokines IL6 in colon and express immunomodulatory effect 75 Inulin is a probiotic that alter the microbiota of the intestine and suppress the inflammation by inhibiting the expression of tumor necrosis factor related cytokines 76 .We employed a multifaceted computational approach to elucidate potential molecular targets for Burdock Inulin, particularly focusing on inflammatory proteins implicated in conditions like DSS-induced colitis.Our findings from molecular docking and simulation shed light on the interaction dynamics between Inulin and key inflammatory proteins including iNOS, COX-2, TNF-alpha, IL-6, and IL-1β.The molecular docking results revealed favorable interactions and binding affinities between Inulin and the selected targets.Notably, IL-1β, COX-2, and iNOS emerged as the most promising targets based on their binding scores and interaction profiles.These findings corroborate previous studies suggesting the potential anti-inflammatory properties of Inulin.

Conclusions
Target validation using in silico computational modeling approaches was employed and validated through molecular simulation-based methods.Our results identified IL-1β, COX2, and iNOS as the best biotargets for inulin.However, simulation-based further validations revealed that among all, iNOS serves as the best target for inulin.The interaction of inulin with iNOS demonstrated stable dynamics, a hydrogen bonding paradigm, and dynamic motion properties.Thus, further experiments are needed to use inulin to target iNOS and reduce DSS-induced colitis and other autoimmune diseases.

Figure 1 .
Figure 1.Interaction paradigm of inulin with iNOS.(a) 3D binding mode of inulin in the active site cavity of iNOS.(b) Show the 2D binding and different interactions in the iNOS-inulin complex.

Figure 2 .
Figure 2. Interaction paradigm of inulin with IL-1B.(a) 3D binding mode of inulin in the active site cavity of IL-1B.(b) Show the 2D binding and different interactions in the IL-1B-inulin complex.

Figure 3 .
Figure 3. Interaction paradigm of inulin with COX-2.(a) 3D binding mode of inulin in the active site cavity of COX-2.(b) Show the 2D binding and different interactions in the COX-2-inulin complex.

Figure 4 .
Figure 4. Interaction paradigm of inulin with TNF-alpha.(a) 3D binding mode of inulin in the active site cavity of TNF-alpha.(b) Show the 2D binding and different interactions in the TNF-alpha-inulin complex.

Figure 5 .
Figure 5. Interaction paradigm of inulin with IL6.(a) 3D binding mode of inulin in the active site cavity of IL6.(b) Show the 2D binding and different interactions in the IL6-inulin complex.

Figure 6 .
Figure 6.Dynamics-based assessment of stability of IL-1b, COX2 and iNOS in bound states with inulin.(a) Show the RMSD results for IL-1b-inulin complex, (b) show the RMSD results for COX2-inulin complex while (c) represent the RMSD graph for iNOS-inulin complex.

Figure 7 .
Figure 7. Residue's flexibility and hydrogen bonding assessment of stability of IL-1β, COX2 and iNOS in bound states with inulin.(a) Show the RMSF results for all the complexes, while (b) show the H-bonds results for all the complexes.

Figure 8 .
Figure 8. Dynamics-based assessment of structural packing of IL-1b, COX2 and iNOS in bound states with inulin.(a) Show the Rg results for IL-1b-inulin complex, (b) show the Rg results for COX2-inulin complex while (c) represent the Rg graph for iNOS-inulin complex.

Figure 9 .
Figure 9. Dynamics-based assessment of trajectories clustering and Free energy landscape of IL-1b, COX2 and iNOS in bound states with inulin.(a-c) Show the PCA results for IL-1β, COX2 and iNOS in bound states with inulin complexes, while (b) show the FEL results for IL-1β, COX2 and iNOS in bound states with inulin complexes.

Table 1 .
Binding free energy results for the three complexes estimated using MM/GBSA approach and reported in kcal/mol.